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Abstract 



We present a phenomenological study of three models with different ef- 
fective degrees of freedom: a Goldstone Boson Exchange (GBE) model which 
is based on quark-meson couplings, the quark derealization, color screening 
model (QDCSM) which is based on quark-gluon couplings with delocalized 
quark wavefunctions, and the Fujiwara-Nijmegen (FN) mixed model which in- 
cludes both quark-meson and quark-gluon couplings. We find that for roughly 
two-thirds of 64 states consisting of pairs of octet and decuplet baryons, the 
three models predict similar effective baryon-baryon interactions. This sug- 
gests that the three very different models, based on different effective degrees 
of freedom, are nonetheless all compatible with respect to baryon spectra and 
baryon-baryon interactions. We also discuss the differences between the three 
models and their separate characteristics. 
PACS numbers: 12.39.-x, 14.20.Pt, 13.75.Cs, 13.75.Ev 







I. INTRODUCTION 



Hadronic interactions are central to strong interaction physics. However, for the time be- 
ing, the fundamental strong interaction theory, Quantum ChromoDynamics (QCD) remains 
too complicated for these to be directly calculated from it. Nevertheless, initial lattice QCD 
calculations of hadronic interactions [I[ have been initiated. Due to the complication of a 
multiquark system, however, quantitative results remain in the distant future. 

Meson exchange models, based on meson-baryon couplings, were developed long before 
QCD |Q and are still the best at fitting the vast collection of NN experimental data |3J. 
Unfortunately, the QCD basis for such effective degrees of freedom in t-channel exchanges 
is not clear at present. Chiral perturbation effective field theory @ is well-based in QCD; 
for a recent derivation, see ||. However, due to the fact that the quark and gluon degrees 
of freedom have been integrated out, it can not be used to study the quark-gluon internal 
structure of hadrons. Nor is it suitable for the study of genuine multiquark-gluon systems 
such as glueballs, hybrids, qqqq systems, and dibaryons. 

R.T.Cahill et al. || have developed an effective field theory which takes spontaneous chi- 
ral symmetry breaking into account. Constituent quarks and Goldstone bosons appear here 
as the effective degrees of freedom for low energy QCD physics. This model has been applied 
both to pion and a meson internal structure and to meson interactions, but not yet to NN 
interactions @. L.Ya.Glozman, D.O.Riska and G.E.Brown || propose a phenomenological 
model, with constituent quarks and Goldstone bosons as the effective degrees of freedom for 
describing baryon spectroscopy and baryon interactions. However, their quark meson cou- 
pling is the linear Yukawa coupling, which reflects neither QCD based chiral perturbation 
nor the nonlinear coupling obtained from the effective field theory referred to above. 

A.Manohar and H.Georgi || have argued that below the chiral symmetry breaking 
scale(~lGev) and within the confinement regime (~lfm), the proper effective degrees of 
freedom are Goldstone bosons, constituent quarks and gluons.(We have constructed an ap- 
proximate QCD derivation of this model, similar to Cahill's approach, but the quantum 
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fluctuations of the gluon field are suppressed.) Such a hybrid quark-gluon-meson exchange 
model has been widely used to describe nucleon-baryon interactions and a semi- quantitative 
fit has been obtained JED] . 

Models with constituent quarks and effective one gluon exchange |TTJ describe hadron 
spectroscopy quite well. The Hamiltonian approach to QCD [|T^] seeks the QCD basis for 
this model. The Schwinger-Dyson equation approach also keeps constituent quarks and 
gluons as the effective degrees of freedom [|l3j. However, direct extension of this model to 
the NN interaction obtains only a repulsive core [ l4| . A variant of this model obtains the 



NN repulsive core and intermediate range attraction simultaneously, by taking into account 
both quark derealization and color screening. This also describes a qualitative similarity 
between nuclear and molecular forces. A semiquantitative fit to the existing NN, NA, and 



NS scattering data and to properties of the deuteron has been obtained as well [15 



The MIT bag model uses current quarks and gluons to describe hadron internal structure 



||16|| . A quark-meson coupling has to be introduced to restore chiral symmetry JTTj to this 
model; Cahill et al. have given an approximate QCD derivation |J. The model has been 
extended to the description of hadron interactions using the R-matrix method and compound 
quark model approach JTj]. 

The Skyrmion model can be viewed as originating from the large N c limit. It has been 
used both for nucleon internal structure and for hadronic interactions The nontopo- 



logical soliton model also acquired an approximate QCD derivation a few years ago f2"0 |. 

There might be even more QCD models that can be listed. But this short review suggests 
that neither theoretical studies of QCD nor phenomenological model analyses can yet deter- 
mine which effective degrees of freedom are best, nor which model is best for the description 
of low energy hadron physics. Therefore, theoretical study of QCD and phenomenological 
analysis are both still needed to extend our understanding of low energy strong interaction 
physics. 

This paper reports a phenomenological study of the effective baryon-baryon interactions 
of the ground state octet and decuplet baryons using three constituent quark models with 



quark-meson, quark-meson-gluon, and quark-gluon effective degrees of freedom respectively. 
The results show that in many channels, (about 2/3 of the total states) the three different 
models give more or less the same predictions. We take this as a phenomenological verifica- 
tion that, even for baryon-baryon interactions, quark models with different effective degrees 
of freedom are quite compatible and that meson exchange effects are modeled, at least to 
some extent, by quark derealization and color screening. In some channels, however, espe- 
cially promising dibaryon states, the different models give drastically different predictions. 
These cannot be viewed as being as reliable as those above; i.e., these results are sensitive 
to model details. Experimental searches in those channels are needed to distinguish the 
validity of the different models. 

This paper is organized as follows: In Sec. II, we describe the Hamiltonians and param- 
eters of the three models. In Sec. Ill, we present the results and discuss them. The final 
section gives our conclusions. 

II. HAMILTONIAN AND PARAMETERS OF THREE MODELS 

From among the models discussed in the introduction, we choose three constituent quark 
models on which to do a comparative phenomenological study; i.e., we calculate the adia- 
batic effective baryon-baryon interactions and dibaryon candidates systematically. The three 
models are: the Glozman-Riska-Brown model, based on constituent quarks and exchange 
of Goldstone bosons (GBE); the Fujiwara model, based on constituent quarks and gluons, 
with the Nijmegen description of one boson exchange (FN); the quark derealization and 
color screening model (QDCSM), based on constituent quarks and gluons but with quark 
derealization and color screening effects included. We choose to consider these three models 
for several reasons: There has been considerable debate recently regarding these effective 
degrees of freedom |21[ ; These models have been extensively studied with respect to baryon 
spectroscopy and baryon-baryon interactions; They can be calculated straightforwardly in a 
consistent, systematic way, using fractional parentage expansions and variational methods. 
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A. Quark and meson effective degrees of freedom: GBE models 

We take the Goldstone Boson Exchange (GBE) model of Glozman's group as the example 
with constituent quarks and Goldstone bosons as the effective degrees of freedom. This 
model has been successfully applied to baryon spectra and also extended to study NN 
interactions. Although recent research shows that scalar and vector mesons should also be 
included to describe baryon-baryon interactions, here we use only the original version, i.e., 
only pseudo-scalar mesons (as Goldstone bosons), which gives a surprisingly good description 
of baryonic spectra. 

In the nonrelativistic meson-exchange quark-quark interaction of any such model, there 
is a 5-function term, which is true only for point-like particles. For particles with structure, 
this 5-function must be smeared out. In the GBE model, two types of smearing were used: 
Gaussian and Yukawa. We designate them as model GBE(a) and GBE(b), respectively. 



1. GBE with Gaussian smearing: GBE (a) 



In this case, the 5-f unction is smeared by a Gaussian function [22]. This version provides a 
successful description of the spectroscopy of baryons. In addition, the short-range repulsion 
of the NN interaction is also obtained. This Gaussian smearing involves two parameters 
corresponding to the short range cutoff point (r ) and the width (a) of a bell-shaped curve. 
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fi n = 139MeV, fi k = 495Me1/, fi v = 5A7MeV, ^ = 958MeV 

m u d = 340M eV, m s = AAOMeV, h± = hl = g ^l = .67 

An An An 

4 s - = 1-206, r = 0.43/m, a = 2.91/m" 1 , V = 0, C = 0.474/m" 2 , b = 0.437/m. 

47T 

Here b is the baryon size parameter (a parameter of the quark orbital Gaussian wave func- 
tion). 



2. GBE with Yukawa smearing: GBE(b) 

Because Gaussian smearing does not meet the well known requirement that the volume 
integral of the pseudoscalar meson exchange interaction should vanish, another smearing 
version, Yukawa smearing, was also applied in the GBE model p3 . 
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m u4 = 3A0MeV,m s = 507MeV,C = 0.77 fm 
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b = 0.537 fm,V = -68Q.AMeV 

The values of b and Vq are obtained by reproducing the iV — A mass difference and the 
nucleon mass with a Gaussian single quark orbital wave function, as is used in standard 



quark model calculations. We also assume a single quark- meson coupling constant, f 2 -, for 



all octet mesons (n, K, if) as the originators of the model have done in their calculations. 
The strange quark mass, m s , is determined by an overall fit to the masses of the strange 
baryons. Since "irrespective of the parametrization, the flavor-spin symmetry is essential in 



the model" 23 1, we presume that we have maintained the essence of the GBE. In addition, 



only the spin-spin component of the pseudoscalar exchange interaction has been included in 
the calculation. We note that there are extensions of the GBE which include pseudoscalar, 



vector and scalar meson exchanges with all possible force components P4" ]. This is expected 
to be necessary for baryon-baryon interactions, but this addition will be left for future 
studies. 



B. Quark-meson-gluon effective degrees of freedom: FN model 

The naive quark-gluon exchange model is quite successful in baryon spectroscopy, but its 
extension to baryon-baryon interactions is much less so. Only a repulsive core is obtained, 
but no intermediate range attraction develops in the best studied NN channel. Meson ex- 
change interactions are introduced into the model to provide this well established attraction 
and the resulting hybrid model with quark, meson and gluon effective degrees of freedom is 
widely used in baryon interaction calculations. There are several different model variations 
in this approach. We choose Fujiwara's model (FN) |25| as an example since it gives a 
simultaneous description of both NN and NY (nucleon-hyperon) interactions. 
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w h = (lfi S cos9 s + fisin9 s \*) x (^cosOs + f§ 'sin8 s \% 
w ff = {-\fisin0 s + ficos9 s \*) x (~f?sirtf s + ficosOgX 1 }) 



m n + = UOMeV, m n = 138M eV, m K = A9QMeV, m v = 547Me V, m v > = 958MeV, 
m 5 = 970MeV,m K = lU5MeV,m e = 80QMeV,m s * = 1250MeV,c 5 = 0.381 
ff = 2.89138, / 8 5 = 1.07509, f[ s = 0.21426, / 8 P5 = 0.26994, 6 PS = -23°, 9 S = 27.78° 
b = 0.616/m, a c = 95MMeV-fm, m ud = 360MeV, a s = 2.1742, A = m s /m ud = 1.526 

C. Quark-gluon effective degrees of freedom: QDCSM 

For this case, we take the quark delocalization, color screening model (QDCSM) as the 
example [IS]. As noted above, the naive quark model with quark and gluon effective degrees 



of freedom gives a good description of individual baryon properties. However, its extension 
to baryon-baryon interactions is not successful. The QDCSM maintains quarks and gluons 
as the effective degrees of freedom, but enlarges the Hilbert space by incorporating quark 
delocalization. It also distinguishes the confinement potential between quark pairs inside 
one baryon and from different baryons by introducing color screening. This is based on the 
recognition that the interaction between quarks from different color singlet baryons must be 
screened as color neutrality is observed to be maintained on a fm scale. To be precise, the 
backflow of color necessary is not explicitly modelled, but is accomplished by an effective 
matrix elements method. The QDCSM gives a reasonable fit to NN, NA and NE scattering 
data. 
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where 4>l{4>r) is the left- (right-) centered quark orbital wave function; the mixing parameter 
e(R s ) is determined variationally for every R s by the dynamics of the six-quark system. 

m u = m d = 313MeV,m a = 634MeV, 

b = 0.603/m,a s = 1.54, a c = 25.UMeV fm~ 2 , p = LO/m" 2 . 

Among these three models, the QDCSM has the fewest parameters. 

In order to do a systematic calculation for all three-flavor (u, d, s) six-quark systems, the 
adiabatic approximation is used to obtain the effective potential between two baryons as a 
function of their separation. The mass of the six-quark system is then estimated as the sum 
of the two separate baryon masses, plus the value of the effective potential at its minimum 
and the zero-point oscillation energy. In the calculation of the six-body matrix elements, 
the fractional parentage expansion is used. The details of this method can be found in Refs. 



T5| , |26|| . Such an adiabatic calculation is reasonable for a systematic survey, but only gives a 
rough estimation of the mass of the six-quark system. A dynamical calculation is needed to 
obtain quantitative results, but is not expected to lead to significant differences at the level 
of accuracy of interest here. 



III. RESULTS AND DISCUSSION 



The baryon-baryon scattering properties calculated are listed in Tables I to VII. Figures 
1-12 show the baryon-baryon effective potentials (in units of MeV) as a function of separation 
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R(fm). 



A. Two smearing versions of the GBE model 

In a series of papers, Glozman et al. used two different versions of smearing of Goldstone 
boson exchanges to calculate baryon masses and the NN interaction. We have extended 
their studies to all possible baryon-baryon S-wave channels consisting of baryon octet and 
decuplet states and find: 

(1) The short-range repulsion obtained from Gaussian smearing is much larger than from 
Yukawa smearing for some states. Taking the (001) deuteron channel as an example, we find 



the short-range repulsion is 2436 MeV for Gaussian smearing ||22|| , but only about 790MeV 
for Yukawa smearing. 

(2) Channel coupling has obvious effects for model(a), but negligible effect for model(b). 
Because of this, some states (about 1/6 of the total) have the similar effective interactions 
in models (a) and (b) for single channel results, but quite different effective interactions for 
results from channel coupling. Such states include (-1,3/2,1) (-1,3/2,2) (-201) (-200) (-213) 
(-221) (-3,1/2,2) (-3,3/2,3) (-402) (-412) , typical examples of which are shown in Fig. 1. 

(3) We identify three groups of states, depending on the minimum value of the potential, 
\Vo\: 

For the first group, the difference of | Vol between models (a) and (b) is about 10 MeV. 
These states include most of the states with small angular momentum (J = 0, 1). These 
states, including (021), (030), (-1,1/2,2), (-1,3/2,0), (-1,5/2,0), (-2,0,2), (-2,2,0), (-5,1/2,0), 
(-6,0,2) and so on, are not sensitive to the different forms of smearing of the ^-function, 
typical examples of which are shown in Fig. 2. 

For the second group, the difference of |Vo| values between model(a) and (b) is greater 
than 20 MeV, but not more than 100 MeV. For these states, intermediate-range attraction 
is always present, but the strength of the attraction varies with the form of smearing. Many 
states with large angular momentum (J = 3) belong to this group, such as (003), (013), 
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(-203), (-3,1/2,0), (-3,1/2,3), (-403) and so on. (See Fig. 2). As for the first group, the 
results do not depend sensitively on the form of the choice of smearing. 

The difference between models (a) and (b) shows up most strongly for the third group 
of states. Comparing the results for (023), (022), (033), (-1,3/2,3), (-1,5/2,3), (-1,5/2,2), 
(-213), (-3,3/2,0), (-3,3/2,3), (-5,1/2,2), (-5,1/2,3), and (-603), (see Fig. 3,) we find a very 
small value for | Vo|, or even no attraction at all, in model (b) , while we find a very large value 
for | Vb|, in the range of 200 to 400 MeV at small equilibrium separation (R s = 0.5 — 0.8/m) 
in model (a). 

In the following comparisons, when we refer to the GBE model, we will mean model (b), 
because it is consistent with the requirement that the volume integral of the pseudoscalar 
meson exchange interaction should vanish. 

B. Comparison among the three different models 

All possible S-wave states of (SIJ) obtained from the ground state octet and decuplet 
baryons in the u, d, and s three-flavor world have been calculated with the three models 
(GBE(a) and (b), FN and QDCSM) based on three sets of effective degrees of freedom. We 
will discuss the general features first and then the details for some interesting aspects. 

1. General features 

(1) Deuteron properties: Since the deuteron is a well known, stable, two-nucleon state 
(SIJ = 001), we take it as a measure of the precision of our model calculation. The three 
different models all show a weak attraction; the minimum of the effective potential is - 
18 MeV (FN), -21 MeV (QDCSM), -3 MeV (GBE), and the separation corresponding to 
the minimum energies is 1.2fm, 1.3fm, 2.0fm, respectively. These reproduce the deuteron 
channel NN interaction qualitatively, but the estimated masses in the three models differ 
by 15-27 MeV from the deuteron energy. The main causes of this deviation are twofold: 
First, the S-D wave mixing due to the one pion exchange tensor interaction has not been 
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taken into account; second the adiabatic approximation and zero-point harmonic oscillation 
energy used in this calculation are rough approximations. A dynamical calculation with 
the extended QDCSM, where the one pion exchange tensor interaction is included, does 
reproduce the deuteron energy, size and D-wave mixing accurately This emphasizes 

that the results reported here only produce a qualitative effective baryon-baryon interaction. 
To study dibaryon states, a careful dynamical calculation is needed. In fact, if the deuteron 
would not have been found by experiment, it would have been difficult for any of the models 
to predict its existence. 

(2) Quark derealization effect: For Fujiwara's quark-meson-gluon coupling model, en- 
larging the model space will not change the results, i.e., systems always prefer to remain 
with small delocalization(e < 0.1) Q, even for those states with minimum energy at small sep- 
arations. A similar conclusion follows for both of Glozman's quark-meson coupling models, 
(a) and (b). On the other hand, for the QDCSM, the quark derealization effect is obvious. 
For some states the derealization parameter e can reach 1.0 at the minimum energy of the 
system. Quark derealization compensates for the explicit meson exchange effect absent in 
the QDCSM. 

(3) Channel coupling effect: In all of these three models, the channel coupling effect is 
small for most cases; the largest effect is about 50-100MeV. For those channels where the 
binding energy is close to the channel threshold, channel coupling may well play a critical 
role in forming dibaryon states. Jaffe's H particle case is a clear example. In those cases, 
the channel coupling should be studied dynamically. 

2. Details 

From the calculated results of three different models, we found that about two-thirds 
of the states within the u, d,and s three-flavor world have quite similar, and in some cases 



1 In Tables III-VII, e is omitted for the FN model. 
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indistinguishable, effective baryon-baryon interactions: 

(1) States with the same, purely repulsive, interaction: 

States with almost the same effective repulsive baryon-baryon interactions include (000), 
(Oil), (023), (022), (033), (032), (-1,3/2,2), (-1,3/2,1), (-1,5/2,3), (-1,5/2,2), (-221), (-403), 
(-411), (-5,1/2,3), (-5,1/2,2), (-603) and (-602). These states have purely repulsive potentials 
in all three models despite their different effective degrees of freedom. There are altogether 
17 states, among all 64 states of 6-quark system, belonging to this case. To simplify the 
presentation, only the four states which have different strangeness are shown in Fig.4. 

(2) States with similar weak attraction : 

(i) Some states have very similar baryon-baryon interactions for all three models, all 
producing a weak attraction, with the minimum energies of the systems appearing at large 
separation. (The separation is a bit smaller for a few states in the FN model). We find 
13 states belonging to this case, such as (001), (010), (021), (-1,1/2,1), (-1,3/2,0), (-201), 
(-211), (-210), (-222), (-3,3/2,2), (-3,3/2,1), (-402) and (-400). Typical ones are shown in 
Fig. 5. Because the binding energies of these states are small, their masses are close to 
the corresponding thresholds. Whether the binding energy is positive or negative is very 
sensitive to the details of the models. It is difficult to determine whether these are strong 
interaction stable states or not, since the precision of the adiabatic approximation is limited. 
A dynamical calculation is especially needed in these cases in order to obtain more reliable 
results. 

(ii) There are states for which the QDCSM and FN give quite similar results, but the 
GBE model yields an even weaker attraction. States belonging to this case include (031), 
(-3,1/2,1), (-401) and (-412), two of which are shown in Fig.6. The QDCSM and FN model 
obtain a minimum of the effective attraction of several tens of MeV, or at most a hundred 
MeV, but the GBE model yields only several MeV for these states. 

(iii) There are states for which the QDCSM and GBE models obtain a weak attrac- 
tion, while the FN model yields a little stronger attraction. These states include: (-1,1/2,0), 
(-1,5/2,1), (-220) and (-3,1/2,0), as shown in Fig. 7. For these states, the QDCSM and GBE 
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models obtain a minimum of several MeV, while in the FN model, the minimum is several 
tens of MeV. 

(iv) For some states, the QDCSM obtains a somewhat stronger attraction than either 
of the other two models. Such states include: (012), (-223) and (-413); the latter two are 
shown in Fig. 8. 

(v) Finally, there are states for which the FN and GBE models obtain slightly dif- 
ferent effective potentials with the results of the QDCSM in between. For the states (020), 
(030) and (-1,5/2,0), the latter two of which are shown in Fig.9, the GBE model yields a 
strong repulsive core and weak intermediate-range attraction; the FN model obtains a weak 
repulsive core and considerable intermediate-range attraction, and the QDCSM obtains, on 
the one hand, a repulsive core similar to that of the GBE model, and on the other, an 
intermediate-range attraction similar to that of the FN model. 

Altogether, we find that the three different models based on the three different effective 
degrees of freedom not only all give a good description of baryon spectra, but they also 
give similar effective baryon-baryon interactions for two thirds of all of the states possible 
within the three-flavor world. This implies that, with regard to baryon spectra and baryon- 
baryon interactions, the models based on the different effective degree of freedoms are still 
compatible. We infer that quark derealization and color screening, working together in 
the QDCSM, must provide a good representation of meson exchange effects at some level of 
accuracy. All of this appears consistent with old conclusions developed from nonperturbative 
QCD: Models with different effective degrees of freedom but with physically equivalent 
results can be obtained for low energy QCD under different approximations. 

Each model has distinct characteristics with respect to the remaining one third of the 
states: 

(1) The QDCSM predicts strong intermediate-range attraction for states with large an- 
gular momentum. As the strangeness becomes more negative, the minimum of the effec- 
tive potential, Vo, becomes less negative, as can be seen from the sequence: (003) with 
V = -359MeV, (-1,1/2,3) with V = -308MeV, (-202) with V = -184MeV, and (-403) 
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where |Vo| is only several MeV. (See Fig.10.) 

(2) The FN model predicts a weak attraction in the effective potential for states with 
small strangeness. Conversely to the QDCSM, here, as the strangeness becomes more neg- 
ative, | Vb | also becomes more negative. For example, for the state (-3,3/2,0) Vo = — 87MeV, 
for (-410) V = -95MeV, for (-5,1/2,1) V = -136MeV, and for (-600) V = -281MeV. (See 
Fig.ll.) 

The Darwin term of the Breit-Fermi interaction has been included in the FN model. We 
find its contribution to the baryon-baryon effective interaction to be small, typically giving 
rise to a weak attraction with a minimum <10MeV. 

(3) Since the GBE model that we used in our calculations does not take into account 
scalar and vector meson contributions, the effective interaction obtained is usually less at- 
tractive than the other two models. One would expect that the GBE will give more similar 
results if these exchanges are incorporated. 



3. Particularly interesting states 

(1) d* dibaryon state (003); (see Fig.10): The QDCSM predicts this to be a tightly bound, 
6-quark state, but the other models predict only a weak intermediate-range attraction, which 
does not provide enough attraction to form a bound state relative to NNtttt threshold. 
Recently, a dynamical calculation in the QDCSM obtained a slightly smaller binding energy 
for the d* (~ 2.18GeV) ||15|| . The difference is due to the fact that the zero-point oscillation 
energy in this case underestimates the relative kinetic energy of the two baryons. However 
the d* mass is still much lower than found in the other models. 

(2) Di-f2 state (-600); (see Fig.ll): The FN model predicts a strong intermediate-range 
attraction with binding energy B = — 62MeV. The mass of the state is 3282 MeV, which is 
lower than the QQ threshold 3345MeV^, so that a strong interaction stable state is predicted. 



Another hybrid model [27[] obtains a similar result; (B is about -100 MeV). The GBE model 
predicts B = — 11 MeV. This mass is slightly lower than the threshold, and is expected 
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to decrease further if scalar and vector meson exchanges are incorporated. The QDCSM 
predicts a mass of 3350 MeV, but a dynamical calculation obtains a smaller mass, which is 
closer to the results of the hybrid models. In this case, the zero-point harmonic oscillation 
energy overestimates the relative kinetic energy. 

(3) H particle (-200); (see Fig. 12): The QDCSM predicts that this is a strong interaction 
stable state with small binding energy, B = —14 MeV; the mass is 2218 MeV. The FN 
model obtains a mass of 2312 MeV, which is higher than the threshold (2231MeV). Zhang's 
dynamical calculation obtains a mass for the H close to the AA threshold p8| . The GBE(a) 
model obtains 2207 MeV, which is a little lower than the threshold, while GBE(b) obtains 
2249 MeV, which is a little higher than the threshold. Almost all of these model results 
are close to the AA threshold. Taken together, they imply that the mass of the H particle 
is sensitive to model details. Experimental searches should (and do) take into account the 
two possibilities of a strong interaction stable state with long lifetime, as well as a narrow 
resonances with a short lifetime. 



IV. CONCLUSION 

The three constituent quark models based on quark-meson, quark-gluon and quark- 
meson-gluon effective degrees of freedom are widely used in the study of hadron spectroscopy 
and hadron interaction, and there is a hot debate about which degrees of freedom are the 



proper ones pl| . For the 64 two baryon channels consisting of octet and decuplet baryons, 
we find that the three models give similar effective baryon-baryon interactions for about 2/3 
of these channels. The three models appear to be different, but in fact are not very different 
in describing effective baryon-baryon interactions. 

For some time, a number of authors have viewed the QDCSM as an exceptional model. 
Some have even claimed that it violates principles of quantum mechanics |[29| . This is in- 
correct, as it is a conventional effective matrix element approach, which we have proposed 
because of our concern regarding a sound theoretical basis for direct extension of the two 
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body confinement interaction from single to multi-hadron systems, where nonlinearity may 
be expected to be significant. Quark model studies of multi-hadron systems provide impor- 
tant checks of the model descriptions of confinement. The results of this paper demonstrate 
that predictions of the QDCSM are not so exceptional as might appear from its definition. 
On the contrary, it embodies meson exchange effects in its own way, i.e., through quark de- 
localization and a different parametrization of quark confinement. Since the effective matrix 
element method is widely used in the Heisenberg version of matrix mechanics, we are at a 



loss to explain how the authors of Ref. |29| have so misunderstood our model as to arrive at 
such incorrect conclusions regarding it. 

This paper also demonstrates that there are differences among these three models. This 



has been shown in baryon spectroscopy ||21|| . It also appears in the predictions for different 
dibaryon masses, especially for states of high spin and high strangeness. More precise hadron 
spectroscopy and experimental study of the dibaryon states predicted by these models are 
needed to distinguish which degree of freedom best describes low energy hadron physics. 

This research is partly supported by the NSF, SED and SSTD of China, and partly by 
the US Department of Energy under contract W-7405-ENG-36. 
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Table I: An asymptotic two-baryon system with channel quantum numbers and delo- 
calization parameter values (e) in GBE(a) and (b) for S — 0, —1, —2. The notation is as 
follows: sc stands for single channel; cc stands for multichannel coupling; B a is the binding 
energy; E a is the mass and Vo is the potential energy at the equilibrium separation, R s ; 
(SIJ) are the strangeness, isospin and total angular momentum quantum numbers for each 
channel. 
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Table I 


GBE(a) 


GBE(b) 




SIJ 




E a 


V 


B a 


e 


R s 


E a 


V 




e 


R s 


Threshold 


0,0,1 


sc 


1889 


-1 


11 


0.0 


2.3 


1891 


-3 


13 


0.0 


2.0 


1878(7VA0 




cc 


1889 


-1 


11 


0.0 


2.3 


1891 


-3 


13 


0.0 


2.0 


0,0,3 


sc 


2503 


-35 


39 


0.1 


0.7 


2474 


-14 


10 


0.1 


1.4 


2464(AA),2158(iViV7r7r) 


0,1,0 


sc 


1888 


-1 


10 


0.0 


2.4 


1892 


-3 


14 


0.0 


1.9 


1878(iViV) 




cc 


1888 


-1 


10 


0.0 


2.4 


1892 


-3 


14 


0.0 


1.9 


0,2,3 


sc 


2333 


-205 


-131 


0.0 


0.8 


2471 


-1 


7 


0.0 


2.6 


2464(AA),2158(7V7V7T7r) 


-1,1/2,3 


sc 


2621 


-40 


5 


0.1 


1.0 


2619 


-24 


2 


0.1 


1.3 


2617(AS*) 
2335(iVA7r7r) 




cc 


2638 


-49 


21 


0.1 


0.8 


2619 


-24 


2 


0.1 


1.3 


-1,3/2,0 


sc 


2141 


-1 


9 


0.0 


2.4 


2144 


-3 


12 


0.0 


1.9 


2132(iV£) 




cc 


2141 


-1 


9 


0.0 


2.4 


2144 


-3 


12 


0.0 


1.9 


-1,3/2,3 


sc 


2523 


-164 


-94 


0.0 


0.8 


2623 


-1 


6 


0.0 


2.5 


2617(AS*) 
2335(iVA7r7r) 




cc 


2523 


-164 


-94 


0.0 


0.8 


2629 


-10 


13 


0.1 


1.4 


-2,0,0 


sc 


2239 





8 


0.0 


— 


2245 





14 


0.0 


— 


2231 (AA) 




cc 


2207 


-106 


-24 


0.0 


0.8 


2249 


-1 


18 


0.0 


1.7 


-2,0,2 


sc 


2481 


-1 


9 


0.0 


2.3 


2487 


-3 


15 


0.0 


1.7 


2472(NS*) 
2397(iVS7r) 




cc 


2481 


-1 


9 


0.0 


2.3 


2487 


-3 


15 


0.0 


1.7 


-2,0,3 


sc 


2785 


-26 


16 


0.1 


1.0 


2779 


-3 


10 


0.0 


1.8 


2770(S*S*) 


-2,1,3 


sc 


2762 


-56 


-3 


0.1 


0.9 


2775 


-16 


10 


0.1 


1.3 


2765(AS*) 
2690(AStt) 




cc 


2701 


-117 


-64 


0.0 


0.9 


2765 


-36 


-1 


0.1 


1.1 


-2,2,0 


sc 


2393 


-1 


7 


0.0 


2.5 


2397 


-3 


11 


0.0 


1.9 


2386(SS) 




cc 


2393 


-1 


7 


0.0 


2.5 


2400 


-3 


14 


0.0 


1.7 


-2,2,2 


sc 


2559 





9 


0.0 




2560 


-1 


10 


0.0 


2.1 


2550(SA) 




cc 


2262 


-471 


-288 


0.0 


0.5 


2568 


-9 


18 


0.1 


1.3 


-2,2,3 


sc 


2635 


-197 


-131 


0.0 


0.8 


2766 





1 


0.0 




2765(AS*) 
2690(AStt) 




cc 


2491 


-393 


-274 


0.0 


0.6 


2769 


-22 


3 


0.1 


1.3 
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Table II: The same as Table I for states with S = —3, —4, —5, —6. 



Table II 


GBE(a) 


GBE(b) 




SIJ 




7— i 

h, a 


T T 

v 


T ) 

-Da 


e 


K s 


7—1 

-&a 


T T 

v 


7 ) 

-Da 


e 


R s 


Threshold 


-3,1/2,0 


sc 


2440 





7 


0.0 




2449 


—1 


15 


0.0 


1.7 


2434(Ar,) 




cc 


2444 


-50 


10 


0.0 


0.9 


2455 


-3 


22 


0.0 


1.4 


-3,1/2,3 


sc 


2867 


-100 


-51 


0.0 


0.9 


2924 


— 1 


6 


0.0 


2.4 


2918(r, E ) 




cc 


2867 


— 100 


—51 


0.0 


0.9 


2930 


—9 


12 


0.1 


1.4 


— 0, o/z, U 


sc 


9QK1 

ZOOl 


1 07 


— 10U 


n n 
U.U 


1 n 
l.U 


9^9Q 
ZOZO 


— 1 


91 
Zl 


n 1 
U.l 


1.0 


2oll(.L>r,J 




cc 


990.D 

zzyu 


— 010 


991 
— ZZ1 


n n 

U.U 


U.4 


Z004 


1 n 

— 1U 


Zo 


n 1 
u. 1 


1 9 
1 .Z 


— 0, of Z, o 


sc 


Zo ( 4 


8 1 

— ol 


— OU 


n 1 

U.l 


n 
u.y 


901 1 

zyi / 


1 1 
— 11 


1 Q 
10 


n 1 
U.l 


1.0 


ZyU4^AAl 6 J 




cc 


9701 
z 1 u 1 


^1 8 


9fH 
— zuo 


n n 

U.U 


U.U 


98Q9 

zoyz 


— 


1 9 
— 1Z 


n 1 

U.l 


1 n 

l.U 


-4,0,0 


sc 


2637 





1 


0.0 




2654 


—5 


18 


0.0 


1.4 


2636(SS) 




cc 


2637 





1 


0.0 




2654 


—5 


18 


0.0 


1 A 

1.4 


-4,0,3 


sc 


3038 


—7b 


—29 


0.0 


0.9 


3073 


— 1 


6 


0.0 


2.4 


3066(S*S*) 


/I 1 9 
-4, 1, Z 


sc 


9898 


^9 
— oz 


9Q 
— Zo 


n n 

U.U 


1 9 
1 .z 


98^9 
ZoOZ 


n 
u 


1 


n n 

U.U 




2ool(^^ ) 




cc 


ZOO 1 


— OOU 


1 oc; 
— iyo 


n n 

U.U 


U.O 


987Q 
Zo ( 


1 Q 

— 10 


91 
Zl 


n 1 

U.l 


1 1 
1.1 


/I 1 ^ 
-4, 1, d 


sc 


9QQ8 
Zyyo 


1 07 
— 1U l 


— oy 


n n 

U.U 


n q 
u.y 


ouoo 


n 
u 





n n 

U.U 




OUO ( ylj i ij 






z,o 1 u 


z,oo 


-Loi- 


n n 

U.U 


U.U 


^Dfi8 


— 1 9 


1 1 

J- 1 


n 1 

U. -L 


1 ^ 


-0, 1/Z,Z 


sc 


9QD7 
Z»U ( 


1 9^ 
— 1ZO 


s'? 
— 00 


n n 

U.U 


1 n 

l.U 


9QQ1 

zyy 1 


n 
u 


1 


n n 

U.U 




zyyu^izj 




cc 


9778 
Alio 


Q99 
— OZZ 


919 
— Z1Z 


n n 

U.U 


U.O 


9008 

zyyo 


u 


7 


n n 

U.U 




-0, 1/Z,d 


sc 


ir\9K 
ouzo 


— zoo 


1 81 
— lol 


n n 

U.U 


n 7 

U. ( 


^907 

ozu 1 


n 
u 


1 
1 


n n 

U.U 




3205(S*fi) 




cc 


3025 


-256 


-181 


0.0 


0.7 


3213 





7 


0.0 




-6,0,0 


sc 


3350 





5 


0.0 


2.6 


3334 


-40 


-11 


0.1 


1.1 


3345(fifi) 


-6,0,1 


sc 


3350 





5 


0.0 


2.6 


3355 


-11 


10 


0.1 


1.3 


-6,0,3 


sc 


3198 


-218 


-147 


0.0 


0.7 


3346 





1 


0.0 
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Table III: An asymptotic two-baryon system with channel quantum numbers and de- 
localization parameter values (e) in the QDCSM, GBE and FN models for S = 0. The 
notation is the same as Table I. (For the FN model, e is omitted.) 
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Table III 


GBE(b) 


FN 


QDCSM 




SIJ 




E a 


V 


B a 


e 


R s 


E a 


V 


B a 


R s 


E a 


V 


B a 


e 


R s 


Threshold 


0,0,0 


sc 


1879 





1 


0.0 


- 


1879 





1 


- 


1879 





1 


0.0 


- 


1878(iViV) 




cc 


1888 


4 


10 


0.0 


- 


1887 


2 


9 


- 


1885 





7 


0.0 


- 




0,0,1 


sc 


1891 


-3 


13 


0.0 


2.0 


1903 


-18 


25 


1.2 


1885 


-20 


7 


0.1 


1.5 






cc 


1891 


-3 


13 


0.0 


2.0 


1903 


-18 


25 


1.2 


1894 


-21 


16 


0.2 


1.3 




0,0,2 


sc 


2478 


-14 


14 


0.1 


1.3 


2479 


-24 


15 


1.1 


2259 


-233 


-205 


1.0 


1.3 


2464(AA) 


0,0,3 


sc 


2474 


-14 


10 


0.1 


1.4 


2479 


-24 


15 


1.1 


2144 


-359 


-320 


1.0 


1.1 


2158(iViV7r7r) 


0,1,0 


sc 


1892 


-3 


14 


0.0 


1.9 


1906 


-35 


28 


1.0 


1893 


-10 


15 


0.1 


1.6 


1878(7VA0 




cc 


1892 


-3 


14 


0.0 


1.9 


1905 


-35 


27 


1.0 


1892 


-9 


14 


0.1 


1.6 




0,1,1 


sc 


1879 





1 


0.0 


— 


1879 





1 


— 


1879 





1 


0.0 


- 






cc 


1885 





7 


0.0 


- 


1890 





12 


0.0 


1885 





7 


0.0 


- 




0,1,2 


sc 


2182 


-2 


11 


0.0 


2.1 


2185 


-14 


14 


1.4 


2122 


-87 


-49 


0.4 


1.2 


2171(iVA) 




cc 


2181 


-4 


11 


0.0 


1.9 


2185 


-14 


14 


1.4 


2122 


-87 


-49 


0.4 


1.2 


2018(jVMr) 


0,1,3 


sc 


2472 


-4 


8 


0.0 


2.0 


2476 


-3 


12 


1.8 


2304 


-188 


-160 


1.0 


1.3 


2464(AA) 


0,2,0 


sc 


2476 


-12 


12 


0.1 


1.4 


2442 


-80 


-22 


0.9 


2383 


-105 


-81 


0.4 


1.4 




0,2,1 


sc 


2181 


-5 


10 


0.0 


1.9 


2170 


-39 


1 


1.2 


2169 


-24 


-2 


0.1 


1.6 


2171(iVA) 




cc 


2181 


-5 


10 


0.0 


1.9 


2186 


-40 


15 


1.0 


2169 


-24 


-2 


0.1 


1.6 


2018(iViV7r) 


0,2,2 


sc 


2172 





1 


0.0 


— 


2172 





1 


— 


2172 





1 


0.0 


— 






cc 


2180 


-1 


9 


0.0 


— 


2172 





6 


— 


2172 





6 


0.0 


— 




0,2,3 


sc 


2471 


-1 


7 


0.0 


— 


2470 





6 


— 


2442 


-41 


-22 


0.6 


1.6 


2464(AA) 


0,3,0 


sc 


2471 


-11 


7 


0.0 


1.6 


2385 


-153 


-79 


0.8 


2416 


-69 


-48 


0.2 


1.5 




0,3,1 


sc 


2473 


-6 


9 


0.0 


1.8 


2440 


-63 


-24 


1.1 


2451 


-32 


-13 


0.2 


1.6 




0,3,2 


sc 


2472 


-1 


8 


0.0 


2.3 


2474 


-3 


10 


1.9 


2474 


-2 


10 


0.1 


2.0 




0,3,3 


sc 


2465 





1 


0.0 




2465 





1 




2465 





1 


0.0 
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Table IV: The same as Table III for S = -1. 



Table IV 


GBE(b) 


FN 


QDCSM 




QTT 

olJ 




TP 






e 


K s 




T T 

Vo 




K s 




Vo 


TD 


e 


K s 


l nresiiolcl 


1 1 /o n 

-1, 1/2,0 


sc 


2U56 


n 



1 


n n 
0.0 




on Kd 
2U56 


n 



1 




2U56 


n 



1 


n n 
0.0 




oncccc f at \ \ 
Zvoo{l\ A ) 




cc 


on 7n 
20 (V 


—I 


1 

16 


n n 
0.0 


1 o 

1.0 


onnn 
2U9U 


Q £ 
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Table VI: The same as Table HI for S = -3, -6. 
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Table VII: The same as Table HI for S = -4, -5. 
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FIGURE CAPTIONS 

Fig.l Channel coupling effect for GBE(a) and (b) where the state is denoted by the symbols 
S,I,J for spin, isospin and angular momentum. The single channel (sc) case for GBE(a) is 
shown as the solid curve; the channel coupling case (cc) for GBE(a) is shown as the thick 
solid curve; the single channel (sc) case for GBE(b) is shown as the dashed curve; the chan- 
nel coupling case (cc) for GBE(b)is shown as the thick dashed curve. 

Fig. 2 The effect of different smearing forms for GBE models. The solid curve is for the 
GBE(a) model; the dashed curve is for the GBE(b) model. 
Fig. 3 Additional channels as in Fig. 2 

Fig. 4 States with similar pure repulsive effective interactions for the GBE, FN and QDCSM 
models. The solid curve is for the QDCSM; the dashed curve is for the GBE model and the 
dotted curve is for the FN model. 

Fig. 5 The same as Fig.4 for states with similar weak attraction for all three models. 

Fig. 6 The same as Fig.4 for states with similar effective interactions for the QDCSM and 

FN models. 

Fig. 7 The same as Fig.4 for states with similar effective interactions for the QDCSM and 
GBE models. 

Fig. 8 The same as Fig.4 for states with similar effective interactions for the GBE and FN 
models. 

Fig. 9 States with effective interactions in the QDCSM which lie between those of the GBE 
and FN model. The curves for each model are as in Fig.4. 

Fig. 10 States with strong intermediate-range attraction in the QDCSM, for which \V \ de- 
creases with the increasing strangeness. The curves for each models are as in Fig.4. 
Fig. 11 States for which \Vq\ increases with the increasing strangeness for the FN model. The 
curves for the models are as in Fig.l. 

Fig. 12 H particle with single channel and channel coupling in the QDCSM, GBE and FN 
models. The curves for the models are as in Fig.4, and the thick curves are for the channel 
coupling case. 
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